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^ ' Reduced gauge theories are theories in which while gauge fields propagate in a bulk, fermion 

fields are localized on a brane. We study dynamical chiral symmetry breaking on a 2-brane 
and a 1-brane in reduced QED3_,_j, and on a 1-brane in reduced QED2_|_i. Since, unlike higher 
dimensional gauge theories, QED3_|_j and QED24.1 are well defined, their reduced versions can 

m ■ 

^ ' serve as a laboratory for studying dynamics in a higher dimensional brane world. The analysis of 

a^ ■ 

^— I in which the conformal symmetry and its breakdown play a crucial role. Explicit solutions of 

, the SD equations in the near-critical regime are obtained and the character of the corresponding 
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. phase transition is described. 



I. INTRODUCTION 



Dynamics in a brane world has recently attracted considerable interest. In most cases, it has been studied in 
higher dimensional theories (for a recent review, see Ref. Q). The aim of this paper is to consider dynamics on 
a brane not in higher dimensions but in a 3+1 and 2+1 dimensional world. More precisely, we study dynamical 
chiral symmetry breaking in the reduced QED3_|_]^ and QED2_(_]^. The term "reduced" implies here that while 
massless gauge fields propagate in a (3+1 or 2+1 dimensional) bulk, fermion fields are localized on a brane. We 
will consider the cases of a 2-brane and a 1-brane in QED3_)^]^ and a 1-brane in QED2+1. We would like also to 
emphasize that though we use the conventional term "QED" for a U{1) abelian gauge theory with fermions, we 
do not specify the origin of the gauge field: it is not necessary electromagnetic field. 

Motivations for considering this type of models are rather obvious. It is well known that relativistic field 
models can serve as effective theories for the description of long wavelength excitations in condensed matter 
systems 0. The reduced QED describes the situation when while fermions are localized on a plane (say, on a 
Cu-0 plane in a high-Tc superconductor) or on a string (polymer like systems), interactions between them are 
provided by a bulk gauge field. Besides that, reduced QEDg^^ can be relevant for the dynamics of cosmological 
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strings g. At last, reduced QED has been recently considered in higher dimensions for the description of the 
mechanism of (quasi)localization of a photon field on a 3-brane . 

Another, more practical, reason is using the reduced QEDg^j^ and QED2_|_]^ as a laboratory for studying 
dynamical chiral symmetry breaking in the brane world. Unlike higher dimensional gauge theories, QEDg,,.]^ 
and QED2_|_i are well defined. While QED3_|_]^ is renormalizable and therefore well defined in perturbative 
theory, QED2^]^ is superrenormalizable and therefore asymptotically free. 

As we will see, the dynamics of chiral symmetry breaking in reduced QED is quite nontrivial. In particular, 
the conformal symmetry (and its breakdown) plays a crucial role in the dynamics. 

The paper is organized as follows. In Sec. |l| general features of reduced QED are described. In Sec. [II 



we study dynamical chiral symmetry breaking on a 2-brane in reduced QEDg^j^. In Sec. IV chiral symmetry 
breaking on a 1-brane in reduced QEDg^j^ is considered. In particular, we discuss subtleties connected with 
spontaneous breakdown of continuous symmetries on a 1-brane. In Sec. ^ chiral symmetry breaking in the 



reduced QEDj^i^j^ with a 1-brane is studied. In Sec. VI we summarize our results. An analysis of the Schwinger- 



Dyson equation in the reduced QEDj.].]^ with a 1-brane is done in Appendix A. 

II. REDUCED QED: GENERAL FEATURES 

In this section, general features of the reduced QED will be described. The QED(-^_]^-)_|_]^ action in Euclidean 
space reads {X = (xd, xi, xjj-i)) 



S = I d^Xi^F^, + - ^^{daA'^)% (1) 



where ^ is a gauge parameter and J° is a fermion current. We will consider the chiral limit (no fermion bare 
mass term) and, for convenience, consistently omit the kinetic term for fermions in the action, restoring it only 
when it is necessary. Integrating over Aa, we get 

3=1 f d^Xd°Yr{X)&°^\x - Y)J\Y), (2) 
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where 



with K — (fcc, fci, ^D-i)- In reduced QED, with a d-brane, we assume that the fermion current has the 
following form: 

r{X) ==0 for a==d+l,rf + 2, 1, 

r{X) ^ f{xD.x^,...,Xd)5''-'-\x) for a = Z?,l,...,d, (4) 
where x = (2:^+1, ■■■,xu_i). Integrating over x and y in Eq.(^, we obtain the reduced d+1 dimensional action 
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[DAeff = \ I d''+'xd^+'yf {x)D'^^,^^^{x - y)f{y), (5) 
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where 



S[Dd]eff = / d'^'^^X 



jFf,^I{-d^)Ft"' +A^jt' + gauge term 



(7) 



^[^V.(- - -7 '^^^^y^ exp - .)) {S,. - (1 - ^ (6) 

with fi,!/ = D, 1, d (the notation for momenta we use here is self-explanatory). 

As it will be shown in the next sections, after integrating over k momenta, the effective action can be rewritten 
in the following general form: 

1 

where is the Laplacian in d+1 dimensions and I{—d^) is a non-local (i.e. integral) operator. 
The following properties of the action (0) are noticeable: 

a) The interacting term A^j^^ is conformally invariant for all D and d. This point will be important for the 
dynamics of chiral symmetry breaking in reduced QED. 

b) When d = D — 2, the kinetic term in expression (^ is finite. However, when d < D — 2, there are 
ultraviolet divergences in it. The reason for that is simple. Because of a delta function in the fermion current 
(^, integrating over x — [xd+i, xd-i) encounters a classical self energy of a point like particle in 13 — d — 1 
dimensions. It is finite in the one dimensional case [d = D — 2) and divergent otherwise. Therefore, when 
d < D — 2, one should regularize the delta function, i.e. introduce a finite thickness for a d-brane. This is an 
additional source of the breakdown of conformal symmetry. Notice that in the reduced QED3_|_]^ with d = 2, 
the kinetic term is both finite and conformally invariant. 

c) Effective action (0) describes fermion fields and a projection of the gauge field on a brane. Since gauge 
bosons can escape from the brane to the bulk, the unitarity does not fulfill in the brane dynamics. In the next 
section, we will discuss explicit manifestations of this feature of reduced QED. 



III. REDUCED QED3+1 WITH A 2-BRANE 



In this section we will study spontaneous chiral symmetry breaking in the reduced QED3_|_]^ with a 2-brane, 
i.e. with D = A, d ^ 2 and k = k^. Integrating over k^ in expression (^), we obtain the bare gauge field 
propagator of an effective 2-f 1 dimensional theory on a 2-brane: 



D 



(0) 



■ exp {ik{x - y)) 5^^ - (1 - C) 



(8) 



where /i, — 4,1,2 and, for convenience, we made the substitution ^ 2^—1, i.e. (1 — ^) 2(1 — 
Introducing a 2+1 vector field A^{x)^ the effective action (|^) can be rewritten in the following form (compare 
with Eq. (0)): 



'[42]e// 



(9) 
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One should add the kinetic term of fermions on a 2-brane to this action: 

Sun = J d^X^Piird^)^. (10) 

As is well known, there are two (two dimensional) inequivalent representations of the Clifford algebra in 2+1 
dimensions. Following Refs. (^-|^, we will consider four component fermion fields which contain these two 
inequivalent representations. In this case, there exists a fermion mass term preserving parity. If there are Nf 
fermion flavors, the symmetry of the action is U{2Nf) The dynamical generation of a fermion mass leads 

to spontaneous breakdown of this symmetry down to U{Nf) x U{Nf). 

A remarkable feature of the action (||) is that it is conformal invariant. Since the initial QED theory is 
renormalizable, one should expect that this feature plays an important role in the dynamics. Our aim is 
to describe spontaneous chiral symmetry breaking in the theory with this effective action. Since there is no 
dimensional parameters in the action (^, a fermion dynamical mass rrid can be induced only through the 
mechanism of dynamical transmutation. In our case, it means that one should intoducc an ultraviolet cutoff A, 
thus breaking the conformal symmetry. Then the dynamical mass, if it arises at all, will be proportional to A. 
We will be especially interested in the near-critical regime of the dynamics, when rrid ^ A. 

The dynamical chiral symmetry breaking in this model is a highly nontrivial problem, and our strategy 
for solving it (at least approximately) will be to find a framework in which the improved ladder (rainbow) 
approximation would be reliable. We recall that while in the ladder (rainbow) approximation there is only a 
Schwinger-Dyson (SD) equation for the fermion propagator (with both the vertex and the gauge field propagator 
being bare), in the improved ladder (rainbow) approximation there are two SD equations (with a bare vertex), 
both for the fermion and for the gauge field propagators. 

The SD equation for the fermion propagator in Minkowski space in the improved ladder approximation has 
the following form: 

G-l(p) = G(°)"\p)+^ J ^^^^G{q)YDy,^^,Ap-ql (11) 



where G^"-* [p) is the bare fermion propagator and 



^42],. = [g,. - (1 - C(fc'))^ j D{k') (12) 

is the full gauge field propagator for which there is its own SD equation (with a bare vertex in this approxi- 
mation). Here we use a general non-local gauge with ^(fc^) being a function of fc^ (a need for considering such 



gauges will soon become clear: see Eq. (21) below). The bare gauge field propagator is now (compare with Eq. 
(i) 

and the full gauge field propagator is related to the vacuum polarization tensor np^(fc): 

^[42],.(fc) - ^[4°C(fc) + "..(fc), n,.(fc) = [g,^ - ^) n(fc2). (14) 
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The structure of the propagator G(p^) is G(p^) = (A{p'^)p - B{p'^)) \ p = j'^Pfj., and from Eqs. and (|l^ 
we obtain the foUowing equations for A{p'^) and B{p^) in the EucUdcan space (po = ipi)'- 



^(^^)-/ (0 ,2^2(;)%(,2) ^[(P-.n(2 + e((p-.)-)). (16) 

Notice that the function D{k'^) is expressed through the vacuum polarization function n(fc^) as 



D{e) = i . (17) 



The full gauge field propagator (|lj) satisfies its own SD equation and therefore is in principle a complicated 
functional of the functions A{p'^) and B{p^). Fortunately, in the present case the situation can be considerably 
simplified. First of all, as we will see below, one can choose a gauge in which the function A(p^) is identically 
equal to 1, and we will use such a gauge. Second, it will be shown that, in the near-critical regime (when 
nid ^ A), the fermion dynamical mass, defined as — 8(171^), is mainly induced in the kinematic region with 
<C fc^. In that region, fermions can be treated as massless, and, if A{p'^) = 1 , the polarization function is 
given by the one loop expression with the fermion propagators of free, 2+1 dimensional, massless fermions. 

For completeness and convenience, however, we will use the one loop expression for n(fc^) taking free fermions 
with the mass m = rrid. On a 2-brane, i.e. in 2+1 dimensions, it is: 

U{k^) = -^[2md + ; ^arctan- ]. (18) 



Notice that 



47r k 2md 



n(fc^) - ^ (19) 

for k 3> rud, and 
for k md- 

When can the improved ladder approximation be reliable? The simplest case would be of course the dynamics 
with a small coupling constant a — e^/47r (notice that a is a bare coupling constant here). In that case, even 
the ladder (rainbow) approximation would be good enough. Unfortunately, as it will be shown below, for small 
a there is no solution with spontaneous chiral symmetry breaking in the reduced QED3_,_]^ with a 2-brane. 
Therefore one should try something else. 

Our initial observation is that the structure of SD equations ( p^ ) and (|l^) is similar to that in usual, non- 
reduced, QEDg^x. It had been recognized long ago that the 1/A^/ expansion can be useful in that theory 
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Though being very nontrivial, the 1/iV/ expansion is helpful in putting under control of nonperturbative 
dynamics. The crucial point is the selection of a "right" gauge in the leading order in 1/-/V/, in which the 
improved ladder approximation would be reliable | pO[ . In particular, appropriate Ward identities have to be 
satisfied in that gauge. In other gauges, the results can be found by gauge-transforming Green's functions 
from the "right" gauge to those gauges. Such a transformation in general changes the initial improved ladder 
approximation to another one, though the gauge invariant quantities remain of course the same. We will adopt 
this strategy for the present problem and, first of all, check the Ward identity for the vertex. Since in this 
approximation, by definition, the vertex is bare, the function A{p^) has to be equal one. It is known [pl| , that 
for the full photon propagator (|l2|), and in arbitrary d space dimensions, this function is identically equal to 1 
if one uses a non-local (in general) gauge with the following gauge function ^(fc^): 

^(^)=d-^^j—^ r dtf'-^Dit). (21) 



z'iD{z) 

We will see that, in the near-critical regime, the momentum region mostly responsible for the mass generation 
i& k':^ nid- As it follows from Eq. (|l9|), in that region Il{k^) — i-e. the function D(k^) (^7|) is proportional 
to k^^ . For such a function D{k^) and d = 2, one gets ^(fc^) — 2/3 (the so called Nash gauge ^), and the gap 
equation takes the form 

''^^^^''^'7 (^?T^7^F=F' ' = 3.^(1 -H^)- ^''^ 

The validity of the Ward identity is a necessary but not of course sufficient condition for the reliability of the 
improved ladder approximation. The crucial point for that is a justification of the use of a bare vertex. This 
approximation for the vertex can be justified in the leading order of the expansion |^-[ic[|. 



Integrating over angles in Eq. (22), we obtain 

Here the ultraviolet cutoff A was introduced. 

Notice that in the momentum region 3> = _B^(m^), the term B^{q^) in the denominator of the integrand 
of expression (|2^) is irrelevant. The only role of this term is to provide a cutoff in the infrared region. Therefore 
one can drop this term, introducing an explicit infrared cutoff in the integral. Then we obtain the following 
equation [x ^ ,y = q^): 

Bi.)^Xr^^Biy/±pl+'-^]. (24) 



The transition from equation (^3|) to equation ( |24| ) corresponds to the so called bifurcation approximation (or 
method) . For the problem of dynamical symmetry breaking, this method was introduced in Rcf. [p^ and since 
then has been widely used in the literature (for a review see Ref. [^). This method is especially appropriate for 
the near-critical dynamics: the closer the dynamics is to a critical (bifurcation) point, the smaller the dynamical 
mass rrid, and therefore the term B^{q^) in the denominator of the integrand (E3), become. 
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It is easy to check that the integral equation (|2J) is equivalent to the differential equation 

x^B" + ^B' + ^B^0 (25) 
with the following two boundary conditions: 

B'iml) = 0, (26) 
(2a;B' + B)U=A2 -0. (27) 



A solution of Eq. (^5|) which satisfies the infrared boundary condition ( |26| ) is 

3/2 / ^ 

B{x)^ Z". ,, sinhP^log^ + ^ , (28) 
a;i/4sinh5 \\ m\ ) ^ ' 

where lo = \/l — 8A and & = \ log yz^i f^nd here we also used the normalization condition B{ni^) — md- The 
ultraviolet boundary condition ( p7[ ) yields the following equation for the dynamical mass: 

th(^logA + j) = _^. (29) 



Obviously, there is no solution rud ^ A for A < Ac- = 1/8. For supercritical values of A (A > Xcr), Eq.(p9[) takes 
the form 

tg(^ log — + arctgi^) = -J/, (30) 
I rud 



where u = VSA — 1. Therefore for small v the mass is 



2lT 

md=iAexp[ h4]. (31) 



The critical line in the plane {Nj, e^) is given by 

where N„iax = Spontaneous chiral symmetry breaking takes place for e > Ccr, and the value N^ax defines 
the upper limit for the number of fcrmion flavors Nf for which spontaneous chiral symmetry breaking is possible. 
When Nf — > N,nax, the critical value e^^ oo. 

Let us now discuss self-consistency of the assumption that the region of momenta q 3> md is mostly responsible 
for the generation of the mass in the near-critical regime (m^ <ti A) . The point is that in this regime In A/m^ ~ 
27r/j/ is large. On the other hand, the behavior of the integrand on the right hand side of equation ( |2^ ) is smooth 
as — > 0. The smooth behavior of the integrand in the infrared region implies that the region < 9 ^5 "^d 
is too small to generate the large logarithm InA/md- It (and therefore the essential singularity in expression 
(^)) is generated in the large region md <C (7 ^ A. A variation of the kernel in the infrared region can at most 
change the overall coefficient in that expression. This heuristic argument is supported by numerical studies of 
integral equation (p3|). 
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The critical line (|3^) implies that it is a strong coupling dynamics, with > e^^, that provides spontaneous 
chiral symmetry breaking on a brane. Indeed, the lowest e^^ corresponds to Nf — 1 and it is e^^ ~ 4.81, i.e. 
acr = elJ4n ~ 0.38. 

This strong coupling dynamics is provided by (essentially) conformal invariant interactions in the most impor- 
tant region of momenta ^ q <C A. Indeed, up to the irrelevant B^(q^) term in the denominator, the kernel of 
integral equation ( p^ ) transforms as K{p^,q^) K{s^p^,s^q^) = s^'^K{p^,q^) under the scale transformation 
p,q ^ sp, sq. This, together with the transformation of the measure (Pq — > s^(Pq, implies that the interactions 
are indeed (essentially) conformal invariant in that region. [In the integral equation (^4|) , the conformal symme- 
try is broken only by the dimensional boundary parameters A and in the integral.] This reflects the presence 
of long range. Coulomb like, interactions which provide the essential singularity in expression (^l|). 

The critical line ( ^2|) corresponds to the so called conformal phase transition (CPT) introduced in Ref. [m) . 
There are the following characteristic features of the CPT: 

a) Unlike the conventional Ginzburg-Landau (GL) phase transition, a parameter governing the phase transi- 
tion in the CPT is connected with a marginal (i.e. renormalizable) operator (in the GL phase transition, such 
a parameter is connected with a relevant (i.e. superrenormalizable) operator; it is usually a mass term). 

b) Though the CPT is a continuous phase transition, there is an abrupt change of the spectrum of light 
excitations at a critical point (line). This is unlike the GL phase transition where the spectrum is continuous 
at a critical point (line). 

In the present model, the parameter governing the phase transition is the coupling constant e. It is connected 
with the marginal operator j^ A^. The spectrum of the light excitations is discontinuous at the critical line 
(^). Indeed, in the subcritical region, with massless fermions, there is a Coulomb, conformally invariant, 
phase describing interactions of massless fermions and gauge bosons. In the supercritical region, with massive 
fermions, there are a lot of bound states, including 2iV| Nambu-Goldstone bosons corresponding to spontaneous 
breakdown of the U{2Nf) to U{Nf) x U{Nf). Therefore these two criterions of the CPT are indeed realized in 
this model. 

Notice that though QED3^]^ is a renormalizable theory, there are new, nonperturbative, divergences in the 



supercritical phase (see Eq. (31)). These divergences are connected not with introducing a 2-brane of vanishing 
thickness in the model but with the strong coupling dynamics. As is well known, such divergences occur in the 
strong coupling phase of QEDg^j^ in the absence of any brane [p"5|~pT|. They lead to breakdown of the conformal 
symmetry (nonperturbative scale anomaly). 

It is instructive to compare the reduced QEDg^.]^ with a 2-brane with the conventional QEDj^x- The SD 
equations in these two models are similar. The difference is in the form of the gauge field propagator. Instead 
expression (p^), one has [0J|]: 

Dik"^) = -jr, , (33) 

where 63 is the (dimensional) coupling constant in QEDj^i and n(/c^) is the (same) polarization function dl^). 
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The appearance of the term A:^/e|, instead makes quite a difference. On the one hand, it provides a 

dynamical ultraviolet cutoff ~ Cg in the SD equation and, on the other hand, since this term is suppressed in 
the region /c <C Cg, it does not contribute to the fermion dynamical mass. This implies reducing screening of 
Coulomb like interactions as compared to the reduced QEDg^j^ with a 2-brane. Indeed, the dynamical mass in 
QED2+1 is 1-0: 

msd ~ 63 exp (^^~~^ ' i^) 

where j/3 = \/8A3 — 1 with A3 — IG/^ir'^Nf. The parameter coincides with v in Eq. only in the limit 
— > cxD, i.e. in the limit of maximally strong interactions in the reduced QED3^2 with a 2-brane.[] 
Therefore we conclude that there are important similarities and important difhrences between the dynamics in 
QED2+]^ and reduced QED3_|^]^ with a 2-brane. Both dynamics are intimately connected with long range Coulomb 
like interactions. Both dynamics provide a realization of the conformal phase transition. In particular, like in the 
reduced QED3^]^, there is an abrupt change of the spectrum of light excitations at the critical point Nf — 
in QEDj^i ||l9|. On the other hand, since QEDj^i is superrenormalizable (and therefore asymptotically free) 
theory, there is no (nonperturbative) ultraviolet divergence in the dynamical mass. Also, its dynamics is more 
effective in generating a fermion mass in that it corresponds to the dynamics in the reduced QED3^j^ when the 
coupling constant e of the latter goes to 00. 

This point is intimately connected with the violation of the unitarity in the brane theory. Indeed, because 



of the first term in the denominator of expression (17), there is an imaginary part for all time like momenta k 
in the gauge field propagator (p^), independently of the value of the mass m^. This feature reflects the process 
of escaping of a gauge boson from the brane to the bulk. This "instability" of brane gauge bosons leads to an 
effective reduction of interactions on the brane. Only in the limit e ^ 00 the gauge bosons are localized on the 
brane, i.e. become "stable". 



IV. REDUCED QED3+1 WITH A 1-BRANE 

In this section we will consider the dynamics in the reduced QED3^]^ with a 1-brane, i.e. with D = A and 
d—1. As it was pointed out in Sec. |l[ there are (classical) ultraviolet divergences in the theory with a 1-brane 
of vanishing thickness in this case. Because of that, one needs to introduce a finite thickness for the brane, 
which will play a role of a regularization parameter. 

To get the reduction 3-1-1^1-1-1, we perform integration in Eq.(0) with the sources taken as 



^The critical value of Nf is iV}'' = 128/37r2 ~ 4.32 in QEDj+i. Since this result was obtained in the framework of the 
expansion, there may be some concern about its reliability [^. Although it would be too strong to say that this 
issue has been finally resolved, different studies indicate that l/Nf corrections are small for Nf around 4 |p||lo|]. 
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r{X) = for a = 2,3, 

r{X) = f{xi, xi)f{x2)f{x3) for a = 0, 1, (35) 

where the regularization function 

/(x) = y''^exp (— ax^), f{x) —!■ S{x), a — > oo. (36) 

Integrating over X2,X3 and 112,1/3 in Eq.(^, we obtain the reduced 1 + 1 dimensional action (|^) with the bare 
gauge field propagator 



1-2 I l2 ; 1,2 I )l2 I jL2 



(27r) 

where k = {ki, ki),k'^ — k\ + kf. It is 



tPk 

2 exp(i/c(x-y))i?['4i]^,(fc), M, 1^-4,1, (37) 



i?(o)(fc2), (38) 



where 



i^(«)(fc^) = -^exp(^)i?^f-^) (39) 



An 2a \ 2a 
and Ei{~x) is the integral exponential function. 

By introducing a 1 + 1 gauge field A^,^ we obtain an effective 1 + 1 dimensional action: 

where i:)(o)'(_a2) ^ i5W(a;)U=_a2. If there were no need for a regularization, this effective action would be 
conformal invariant. The finite thickness of the 1-brane breaks the conformal symmetry. We will return to this 
point below. 

One should add the kinetic term of fermions on a 1-brane to the action (^): 

Smu = y" d'xi,{i^^d^)^. (41) 

We will consider N j two component (i.e. vector like) fermions. The chiral group is U{J>i j)h x \J(N f)]^. Our 
aim is to find whether a fermion dynamical mass is generated in the theory with effective action (|40|). Naively, 
one might expect that in this case the chiral symmetry U{Nf)]^ x U{Nf)ji would be spontaneously broken 
down to its vector subgroup U{Nf)v- However, this is not the case in 1+1 dimensions. Due to the Mermin- 
Wagner-Coleman (MWC) theorem |2^, there cannot be spontaneous breakdown of continuous symmetries in 
1 + 1 dimensions. The MWC theorem is based on the fact that gapless Nambu-Goldstone bosons cannot exist 
in 1 + 1 dimensions. It however does not prevent a generation of a fermion mass. In this case, the so called 
Berezinski-Kosterlitz-Thouless phase would realize [^l],^ . We will return to this point at the end of this section. 
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As in the previous section, we will change the initial gauge in such a way that the full gauge propagator takes 
the form: 



(42) 



Our aim is to find such a function ^(fc^) that the fermion wave function A(fc^) would be identically equal to 
1 in the improved ladder approximation. Fortunately, as Eq. ( |2l] ) shows, for d = 1 the choice of ^ = 1, i.e. 
the Feynman gauge, provides A(fc^) — 1 for any gauge field propagator. Then the SD equation for the fermion 
mass function takes the following form in the improved ladder approximation: 



cPkBiJP) 



■D{{p-kf). 



(43) 



fc2+B2(fc2) 

As was shown in the previous section, in the reduced QEDg^j^ with a 2-brane, there is no generation of the 
fermion mass for a small coupling constant a, and we used the 1/-/V/ expansion in order to justify the use of 
the improved ladder approximation. As we will see below, the situation in the reduced QED3^]^ with a 1-brane 
is different: there may exist a solution with a nonzero even for an arbitrary small a. For such a a, one can 
expect that even the ladder approximation is justifiable. 

Using the same arguments as in the reduced QED3_|_]^ with a 2-brane, one can show that in the present case 
the momentum region yielding a dominant contribution to SD equation ( ^ ) is ^ fc <C A (we will see 
below that the parameter a}/'^ (the inverse thickness of the 1-brane) plays the role of an ultraviolet cutoff A). 
Therefore, like in the previous section, for the vacuum polarization function one can use the one-loop expression 
with the propagators for free massive fermions with the mass md- In 1 + 1 dimensions, this expression is: 



n(fc2 



N 



f 



1 



2ml 



In 



v/fc2 -I- 4m2 + ' 



VF(fc2 + 4^ v/fc2 + 4m2 



The asymptotics of this expression are: 

n(fc2) 

n(fc2) 

From Eqs. ( 

D{e) 



6? ' 



for k^ ^1712, 



f 



for k^ » mi. 



and ( p^ we find that for momenta k^ ^ 2a the photon propagator (42) rapidly decreases: 

1 2ae2 



2e^a 



^ 47rfc2 ■ 



(44) 

(45) 
(46) 

(47) 



Therefore the whole integrand in SD equation ( [4^ ) rapidly decreases for k^ ^ 2a. Because of that, one can 
neglect the region of those large momenta and put the cutoff A2 = 2a or, more precisely, A2 = 2aexp(— 7), 
with 7 the Eulcr constant, in the SD equation. 

This also implies that one can keep the leading order term in expansion of D^'^^ (|39|) in k'^/a: 



Do{k^) ^ — lo{ 



2a exp (—7) 



An "° fc2 
From here and Eqs. (^), (^), and (^ ) we find that 



e2 A2 
4^^°S#- 



(48) 
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An 



D-\k^)^ at fc2< 



e2 log ^ 



i?-i(fc2)c.-i!^ + ^, at e>ml (49) 
log -j± 

Now we proceed at solving the SD equation. In order to get a hint of the character of the solution, we will 
first consider the so-called a constant mass approximation, taking the external momentum being equal to zero 
and replacing the running mass function in the integrand by its value md — B{0). Then we get the equation 

_ f d^k 1 2 



, (27r)2 fc2 + ml 42_^ _^ Nj_ 

e2 log -J- ^ 

The main contribution comes from the range of momenta fc2 3> rn^. Therefore one can omit the term m2 in 
denominator and put instead the parameter as the lower limit in the integral. Then the integral can be 
easily evaluated and one gets the following algebraic transcendental equation: 



2Nf 



1 f, A2 47r2 , 
l = 777r log-|--^log 



e'^f , A2 



47r2 



(51) 



Introducing the variable y — {e Nf /An ) log(A2/m^), it can be rewritten as 



— ^ =2/-log(l + 2;). (52) 

ZTT 

The function on the right hand side of this equation is monotonically increasing, starting from zero value at 
2/ = and going to cx3 as y — > oo. Therefore this equation always has a solution. However, the character of the 
solution depends on the value of the parameter e^N'j /2tt^. Indeed, one gets: 

ml^kle^v(^-^^, for ^ « 1, (53) 

ml^Kl[-^\ exp(-2iV/), for ^ » 1- (54) 

Notice that solution ( ^ ) corresponds to a weak coupling regime. It does not depend on Nf and therefore comes 
from the range of momenta in the integral equation where one can neglect the vacuum polarization, i.e. one 
can use the ladder (rainbow) approximation in this case. 

A closer look at Eq.(^l]) reveals that the solution ( p3| ) emerges from the range of momenta where a double 
logarithmic contribution dominates in the gap equation (|50|). On the other hand, the solution at large Nj (Eq. 
(^)) comes from the region of momenta generating a one logarithmic contribution. 

We will turn now at studying SD equation (^3|) for the running mass function. In order to integrate over the 
angle variable there, we will use the following conventional approximation for the vector boson propagator 
(for a review see Ref. [Q): 

D{{p ~ kf) ~ D{p^)9{p^ - fc2) + D{k^)e{k^ ~ p^). (55) 

This approximation is justifiable: the measure of the only "dangerous" (for this approximation) region, with 
|p2| ~ |fc2|^ jg gjjja.ll and the dependence of the propagator D{{p — k)'^) on the angular variable is rather smooth. 
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Neglecting then the term B^(k^) in the denominator and instead putting the infrared cutoff = i3^(m^) in 
the integral (the bifurcation approximation discussed in the previous section), one gets a simple integral equation. 
It is easy to show that it is equivalent to a differential equation with two (infrared (IR) and ultraviolet (UV)) 
boundary conditions: 



(56) 



B'{x)\,^^. = 0, [D{x)B'{x) - D'{x)B{x)] 



0. 



It is convenient to introduce the variable z = {e^ N j / in^ ) log(A^/x) in terms of which Eq.(^6|) becomes: 
B"{z)-—^B'{z)-^^D'{z)B{z)=0, D{z) 



D'{z) e^Nf 
Together with the boundary conditions, it can be rewritten as 

27r2 



Nf z + 1' 



B"{z) + ^—B'{z)+ 
B'{z) 



B{z) ^ 0, 



.2 =0, [zB'{z)-B{z)] 



0. 



(57) 

(58) 

(59) 
(60) 



The solution B[z) satisfying the UV boundary condition and the normalization B[x = rnjj) ~ irid is given by 
B{z) = md 



zo + lV^' sinh[f log(z + l)] 



z + iy sinh[f log(zo + l)] ' 



87r2 



where zq = z{x = rn^). The IR boundary condition leads to the equation for the dynamical mass: 



tanh 



2 log(^o + l) 



■ to. 



For real lo {e Nj/Si: > 1) it can be easily solved: 



For large e^Nj this solution becomes 



l+OJ 
1-UJ 



1/u 



- 1 



(61) 



(62) 



(63) 



2 ]\t2 \ c^JV, 



/ 



exp {-2Nf ) , 



(64) 



and coincides, up to minor difference in preexponential factor, with expression ( p^ obtained in the constant 
mass approximation. 

The line a = 6^/471 = 271 /N'j divides the region in plane (a, Nf) in two parts with different dependence of a 
dynamical mass on a and Nf. Indeed, at a < 2Tr/Nf, when uj = iv,v — ./STr^/e^A^? — 1, Eq.(|65) gives 



- K exp 



(65) 
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The ratio of powers of two exponents in ( p^ ) and ( p3| ) is « 1.11 that shows that the constant mass 

approximation in this case is also rather rehable. It is pecuUar that the expression ( |65| ) for a dynamical mass 
coincides with the expression for a dynamical mass generated by a magnetic field in quenched QEDg,,.]^ (see 
Eq. (Ill) in Rcf. |24[ ). In fact, in the ladder (rainbow) approximation, used in the weak coupling regime, the 
present SD equation essentially coincides with the SD equation in that paper (see especially Appendix C there). 
The origin of this similarity is in the dimensional reduction 3 + 1 — > 1 + 1 in the dynamics of spontaneous chiral 
symmetry breaking in a magnetic field ||2^ . 

The existence of the two types solutions, corresponding to the weak (e^ = e^Nf ^ 1) and the strong (e^ 3> 1) 
coupling regimes, is intriguing. While the strong coupling solution essentially coincides with that in the 1 + 1 
dimensional Thirring model (see below), the weak coupling one yields a new type solution, characteristic for 
a 1-brane physics in a 3 + 1 dimensional bulk. These two solutions are generated by very different dynamics: 
while in the strong coupling regime the gauge field propagator is dominated by the 1-brane vacuum polarization 
operator, in the weak coupling one the propagator is dominated by the bare term coming from the bulk. In 
particular, while the polarization operator is generated by the conformal invariant interaction jfj,A^ , the bare 
term breaks the conformal symmetry as result of a finite thickness l/a^^^ of a 1-brane. We will argue below 
that this point can be important in the connection with the MWC theorem. 

As we already stated above, there cannot be spontaneous breakdown of a continuous symmetry in 1 + 1 
dimensions (the MWC theorem) This happens because strong fiuctuations of would be NG modes lead 

to vanishing order parameter connected with such a breakdown. Let us recall how this theorem is realized 
in the case of the 1 + 1 dimensional Thirring model with the color group U{Nc) and the chiral symmetry 
U{Nf)L X U{Nf)ii. It is relevant for our case since the dynamics of the strong coupling solution found above 
essentially coincides with the dynamics of the Thirring model with the color group U{\). Indeed, since for 
this solution the gauge field propagator is dominated by the 1-brane vacuum polarization function, which is 
essentially constant in this case (see Eq. (p6|)), the interaction is of a current x current form, as in the Thirring 
model. 

First of all, using the Fierz identities, it is easy to show that, in 1 + 1 dimensions, the Thirring model is 



equivalent to the Gross-Neveu (GN) model |25|. The interaction term of the latter is: 

5f„f = j d'x^ [i^X^^f + i^XHj'^Pr] , (66) 

where A*' are flavor matrices, s = 0, 1, ■■■iNj — 1, and the summation over s and color indices of the fermion 
fields is assumed. The A'^ matrices are normalized according to tr(A''A'^) = 25'^^. 

Let us first consider the case of the U{1)l x U{1)r chiral group. In this case the model is soluble There 
is a nonzero dynamical mass for fermions for all Nc > 2. However, there is no NG boson in the model. Instead 
of that, there is a Berezinski-Kosterlitz-Thouless (BKT) gapless mode. This mode is described by the exponent 
field U{x) = exp{i6{x)), where 9 satisfies the constraint < d{x) < 2tt. More precisely, the BKT mode is 
described by a usual Lagrangian density of a massless free field, f /2{d^9d^6) with / ~ Nc/Att. However, the 
corresponding observables are described not by Green's functions of the field 6 but by Green's functions of the 
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field U{x) and its derivatives, including the derivative dfj,9 = iUd^Jj'^ . The point is that while the propagator 
and other Green's functions of the 0{x) field do not exist in 1 + 1 dimensions (they are divergent for all x), 
Green's functions of the U {x) field and its derivatives are well defined. Moreover, the corresponding field theory 
is conformal invariant and the parameter / defines anomalous dimensions of its Green's functions. 

The case of the GN model with one color is special. For Nc = \ and the chiral group U{\)l x U{1)r, fermions 
are massless and, moreover, the bosonization of the model leads exactly to the Lagrangian of the free massless 
BKT mode |p7| . Therefore in this particular case, the whole dynamics is conformal invariant. 

Though the dynamics with Nf > 2 is more involved, some of the basic points described above survive. In this 
case one should distinguish the U{\)l x U{\)ii and the SU{Nf)L x SU{Nf)ii sectors. The dynamics in the first 
one is essentially the same as in the model with Nf — 1 and one should expect that while for Nc> 2 fermions 
are massive, they become massless for Nc = 1. In the second sector, because of a strong self- interaction between 
Nj — 1 would be NG bosons, all of them acquire a (same) mass, thus leading to a Wigner realization of the 



dynamics with the exact SU{Nf)L x SU{Nf)f{ symmetry |28 



In our case the number of colors Nc = 1. Docs it necessarily imply that the dynamical mass of fermions will 
disappear in the exact solution in the reduced QED^^i with a 1-brane? We do not think that the situation is 
so simple. First of all, even the status of the Goldstone theorem is not completely clear in this case: some of 
the assumptions the theorem is based on are violated in the brane world. Indeed, in the initial bulk theory, 
the (_D — 1) + 1 Lorentz invariance is broken because of the presence of a d-brane. On the other hand, while 
on a d-brane the d + 1 Lorentz symmetry is preserved, the corresponding effective theory in nonlocal. Second, 
as was emphasized above, in the 1 + 1 dimensional Thirring (or Gross-Neveu) model, it is important that the 
conformal symmetry is exact in the sector with the BKT field U{x). On the other hand, in the reduced QEDg,,.]^ 
with a 1-brane, the conformal symmetry is necessarily broken by a finite thickness of the brane. The latter 
is especially important for the weak coupling solution ( |65|) in which the gauge field propagator is dominated 
by the the bare term ( psf ) which explicitly breaks the conformal symmetry. The dynamics described by that 
solution is very different from that of the Thirring model. 

It remains a challenge to clarify these various issues in the brane dynamics. 



V. REDUCED QED2+1 WITH A 1-BRANE 

In this section we will study spontaneous chiral symmetry breaking in the reduced QED2^x with a 1-brane, i.e. 
with D=3 and d=l. Recall that the gauge coupling constant is dimensional in 2 + 1 dimensions: its dimension 
is [63] — mi, and we will see that the parameter e'^Nf plays the role of an ultraviolet cutoff, which is a typical 
feature for Notice also that as it follows from the discussion in Sec. |^ there is no need for 

introducing a finite thickness for a 1-brane in a 2 + 1 dimensional bulk. 

With trivial modifications, the effective action can be derived as in the case of the reduced QED3_|_;^ with a 
2-brane. It is: 



15 



'5'[32]e// ^ I d X 



(67) 



As in the previous section, wc will consider Nf two component fcrmion fields (see Eq. ([4l|)). The chiral group 
is UiNf)L X U{Nf)R. 

The full photon propagator in a nonlocal gauge is given by 

1 



fc2 



2k 



n(fc2)' 



(68) 



where the vacuum polarization function n(fc ) is given in Eq. ([44|) . As was shown in Sec. IV, the convenient 
choice of the gauge for the study of fermion dynamics on a 1-brane is = 1. Then the SD equation for the 
fermion mass function takes the form of Eq.(^) with the function D(k) given now by 

Like in the case of the reduced QED3_|_]^ with a 1-brane, the bare term 2k /e^ breaks the conformal symmetry. 
However, this bare term is very different from that one in Eq. (|5^). Its strong dependence on momentum 
implies that it is important both in the infrared and ultraviolet regions. Taking into account this term and 
the asymptotics of the vacuum polarization function Il{k^) (^) and (|46|), one concludes that the dominant, 
logarithmic, contribution to the SD equation should come from the range of momenta 127rm^/e§Af/ < k < 
elNf/2TT. 

An analysis of this SD equation is done in Appendix A. It is shown there that a solution with a nonzero 
dynamical mass exists for all values of Nf and Cg. It is also shown that the dynamical mass satisfies the 
following constraint: 

In the case when Nf ^ 1, the dynamical mass is: 

md^^exp[-(7V/ + -+7-31og7)]. (71) 

It is interesting that, unlike the previous model with D = A and d = 1, in this model the constant mass 
approximation is unreliable. In particular, it is not difficult to show (see Appendix A) that it would yield the 
following expression for the dynamical mass in the case Nf 3> 1: 

which is very different from expression (^Tj). The reason for that is that, unlike the previous model, the bare 
term in the propagator ( |69| ) (now strongly depending on momentum) does not decouple even in the dynamical 
regime with Nf ^ 1. 

We would like also to add that all remarks made in Sec. |^ concerning the status of the problem of the 
fermion mass generation on a 1-brane, in particular, its connection with the Mermin-Wagner-Coleman theorem, 
are also relevant for the present case. 



16 



VI. CONCLUSION 



The dynamics of chiral symmetry breaking in reduced QED is rich and quite nontriviaL Its characteristic 
features are intimately connected with the structure of the gauge field propagator. It includes two terms: the 
vacuum polarization function, completely defined by the brane dynamics, and the "bare" term coming form the 
bulk. The vacuum polarization function is connected with the conformal invariant term jf^A^. Therefore, since 
in 1 + 1 and 2+1 dimensions there are no divergences in the polarization function, it is conformal invariant for 
massless fermions. This feature essentially survives in the near-critical regime of chiral symmetry breaking: in 
this regime, a fermion dynamical mass is small and the dominant region is that with momenta k 3> rud- 

On the other hand, in many cases, the bare term breaks the conformal symmetry: either because of a finite 
thickness of a brane or because an initial bulk theory (as QEDj^x) conformal invariant. The interplay 

between those two dynamical sources provides rich nonperturbative dynamics. 

In this paper, the improved rainbow approximation (with a bare vertex) was used. It would be interesting to 
study the dynamics beyond this approximation, though it is not straightforward at all. The point is that, besides 
the bare spin structure 7^, the vertex can have other ones. For example, in the case of a 2-brane, there are in 
principle 11 other structures. The crucial point in the present analysis is decoupling of the Schwinger-Dyson 
equations. Therefore the role of the gauge where the function A{p^) = 1 is very important. As it is discussed in 



Sec. [II, for the bare vertex, one can find such a gauge for any vector boson propagator. However, beyond the 
approximation with the bare vertex, new structures in the vertex can appear. In order to find these structures, 
one either should consider the equation for the vertex (that is quite complicated) or try to construct an ansatz 
for the vertex consistent with such general constraints as Ward identities, the absence of kinematic singularities, 
the correct perturbative limit, etc.. This last approach was successful in 3+1 dimensional QED ]30| ]. However, 
studies of this problem in 2+1 dimensional QED (which is similar to the dynamics on a 2-brane) have revealed 
that it is a hard (and still unresolved) problem We hope to turn to this problem elsewhere. 

At last, we would like to indicate that this analysis can be useful for studying dynamical chiral symmetry 
breaking in higher dimensional brane theories p2| , ^ . In this connection, it is noticeable that in Ref. ||3^ 
the consequences of the existence of a ultraviolet stable fixed point in higher dimensional gauge theories were 
considered. 
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APPENDIX A: ANALYSIS OF THE GAP EQUATION FOR THE REDUCED QEDj+i WITH A 

1-BRANE 

In this Appendix we analyse the SD equation for the case of the QEDg^x with a 1-brane. The equation has 
the form 

where 

D{{p-kf 



and the vacuum polarization function n(fc^) is given in Eq. (p4[). 

We will first obtain the constraint (|70| ) for the dynamical mass. We begin by deriving the lower limit for m^j. 
As was already indicated in Sec. the dominant contribution to SD equation (Al) comes from the range of 



momenta p, < k < A with infrared and ultraviolet cutoffs given by /i = 127rTO^/A^ye§, A — e'^Nf/2TT. Since the 
kernel of this integral equation is positive (corresponding to an attractive interaction), we obviously obtain a 
lower limit for rrid if integrate only over this range of momenta and, furthermore, replace D{{p — fc)^) in the 
kernel by its minimal value in this interval. Taking into account Eq. (^sf), one finds that the minimal value is 
Then the gap equation becomes simple: 

2x 2 d^kB{k^) TT 



It has the following solution: 



Bip^) ^md~ -^-^ exp{-2Nf). (A3) 



Since the initial interaction is stronger, the true is larger than the value ( AS ). 

Let us find an estimate from above for the dynamical mass. To do this, we consider the integral equation at 
p2 = 0. It is 

2 f cPkBik'^) ,90, , , , 

^(°)-(2^yi:^^(^'"^^)' 

and we explicitly indicated the dependence of the interaction kernel D(k^,rn^) on the dynamical mass m^. 



Eq.(A4) is equivalent to 



^ = 77^ i'^P^D{k\m% (A5) 



(27r)2 J P 

where /(fc^) = ^B(ti)- follows from the gap equation that B'{p^) < 0, i.e., B{p^) is a decreasing function of 



Therefore, f{k^) < 1 in Eq.(A5) for k^ > 0. In the case of the constant mass approximation (where B{k^) 



is a constant B{k ) = M) the square mass M satisfies the following gap equation: 
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1 = 



(27r)2 J k^ + M^ 



D{k'^,M^ 



(A6) 



By using asymptotics (|45| ) and ([46| ) for the vacuum polarization function, the integration region in Eq.(A6) is 
divided into two regions k < MVG and k > MVQ R that gives us the following gap equation: 



1 

1 = - 



MV6 



fc2 +7\/2 2 , Wj._fe 
' 67r 



dk 



(A7) 



We can further neglect /c^ term in comparison to in the first integral in Eq.(A7), while in the second one 
we can neglect in comparison to k^. Evaluating the integrals, we come to the following expression: 



1 



6 



■ log 1 , 



1 



Nf4 



log 1 , 



The corresponding solution for a small dynamical mass (Af <C eg) is 

M 



Nfel 



■ exp 



7 



(A8) 



(A9) 



2V67r 

It is obviously valid for A^/ ^ 1. 

Let us prove that m?^ < M"^, where is the solution of the gap equation with the running mass function, 
by assuming the opposite and then showing that it leads to a contradiction. 

So let us assume that > Af^ and consider the integral 

d'^kfiP'' 



fc2 + 



D{k',mi). 



Since /(A:^) < 1, we have 



By calculating 



d^kfik"^) 
k"^ + 



D{k\ml)< 



d^k 



A;2 + m? 



Dik^ml). 



(AlO) 



(All) 



d^k 



fc2 



rD{k^,mj), 



one can show that I'{m2) < 0, i.e., /(m^) is a decreasing function of m^. Since we assumed that > M , we 



have 



and, consequently, we obtain that 



d^k 



fc2. 



-D{e,ml) < 1 



fc2 + 



D{k\mi)<l 



(A12) 



^The value here was determined from matching small and large k asymptotics of the vacuum polarization function. 
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for all > M^. Then, since we cannot satisfy equation (A5) with rn^ > M^, the assumption that > 
leads to a contradiction. Therefore, we get the inequality < with M given in Eq. (A9). This and the 
lower limit we obtained earlier lead us to constraint 



Can one get an explicit solution of the integral equation (Al) in a reliable approximation? The answer is 
affirmative. 

To solve Eq.( |Al[ ), we use approximation ( ^5| ) in order to be able to perform integration over angles. Since we 
already know that the main (logarithmic) contribution comes from the range of momenta 12nm'^/Nfe'^ < k < 
e'^Nf/2TT, we put infrared and ultraviolet cutoffs in the integral equation at /i = 1277771^/ NfC^ and A = e\N f /2tt^ 



respectively. Then the integral equation (Al) takes the form: 



B{p) 



D{p) 



dkkB{k) f dkkB{k)D{k) 



k^ 



k^ 



(A13) 



Furthermore, we approximate the function D{p) on the interval /i < p < A as 

1 



2p , AT/ ' 



(A14) 



where the parameter p,„ is determined from the condition of continuity of D{p) at the point p = pm, {pn 



It is convenient to represent the mass function as 



6m^ « md\/&). 



B{p) = B,{p)0{prn -P)+ Buip)0{p - Pn 



(A15) 



For "infrared" Bi and "ultraviolet" B^ (with respect to the parameter pm) parts of the mass function, we get 
the following equations: 

P Pm A 



6m^ f dkkBi{k) 



Nfp"^ 



dkkB,{k) Qml 1 
k'^ + Nfk"^ TT 



dkkBu{k) 
k^ + rnj, 



1 



2k , Nf_ 



(A16) 



Buip) = - 



2P 1 ^7 fc2. 

^ Pn 



dkkBuik) f dkkBu{k) 



fc2 



Taking the derivatives on the both sides of these equations, we get: 



B'M = 



B'uip) 



l2TTml f dkkBi{k) 



Nfp^ 



fc2 + : 



dkkB^(k) 



1 2 



dkkB^ik) 

fc2 + TO^ 



dkkBuik) 
fc2 + m? 



Differentiating the last equations once more time we obtain 



, 3 , 12771^ 

B'{p) + -Bp + — 2^ 

' ^ ' p ' Nf p^{p^ + ml) 



B{p) 



B':{P) 



2lT 



-B'M 



pBu{p) 



27r 



{p+'-^?{p^+ml) 



= 0, 
= 0. 



(A17) 



(A18) 
(A19) 

(A20) 
(A21) 
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We have also the following IR and UV boundary conditions: 



Blip) =0, [(p + A)i?„(p)]' 

P=A' 



0, 



(A22) 



where the prime denotes the derivative with respect to p. Furthermore, the mass function is continuous at the 
point prn, therefore, Bi{pm) = Bu{pm) and the first derivatives satisfy 



(A23) 



(the condition of continuity and (A23) follow from Eqs.(A16), (A17) and Eqs.(A18), (A19), respectively). 



The general solution of Eq.(A2C) is given in terms of hypergeometric functions: 



B{p') - Ci ^ F 



2 \ — 



-l+UJ l+UJ ^ p^ 
n '^^'1+^; 2 

2 2 ml 



where lo = ^ \ 



12 

Nt 



1 + 1 — 
2 2 TV^ ' 



The IR boundary condition gives a relation between the constants Ci and C2 



(A24) 



l+UJ l+UJ 



2 ' 2 

\ — UJ 1 — UJ 



rrid 



+ C2(1+W) 

= 0, 



where we used the formula for differentiating the hypergeometric function |29| 



[^a+n-i^ (a, b; c; z)] = {a)nz''-^F {a + n, b; c; z) . 



(A25) 



(A26) 



Since for Bu{p) the corresponding momenta are larger than [p > p„i > rrid), we approximate p^ + by p^ 



in Eq.(A21). This gives us: 



P ' 2rr 
2 



4Nf " 



0. 



Introducing the variable z = ~2'Kp/Nfe\ and making the substitution 



z- 1 



Eq. (A27) reduces to the hypergeometric differential equation 



A solution regular at zero is 



(A27) 



(A28) 



(A29) 



(A30) 



and the second independent solution is p9[ 
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/2(t) = (l-t)F(i±^,i-^;2;l-< 



Thus, the general solution for the mass function is 

„ z f 1 + V 1 — V z 
Bu{z) = C^z F { ^— , ^— ; 2 



l-z 



z-1 



1 fl + v l-v 1 
Ciz F { ^— ;2 



l-z 



1 - z 



The UV boundary condition (A22), which can be rewritten as 



0, 



allows us to fix the ratio C3/C4 

P ( 1 + ^ l — u . O. TP ( l+i^ l — iy . 1 . j,"! 

L^S _ V 2 ' 2 ' ' 2/ V 2 ' 2 ' ' 2y 



C4 F (i±il, i-^; 1; -i) + F (1±^, i-^; 2; -1) ' 
where the formula for differentiating the hypergeomctric function 

— [z^~'F (a, b; c; z)] = (c - n)„z=-i-"F {a, b:c-n; z) 



(A31) 



(A32) 



(A33) 



(A34) 



(A35) 



has been used. For iV/ > 1 we have C3/C4 ~ 1/8A^/. 

Finally, matching the solutions Bi [p) and i?„ [p) at the point p = p^ — m\/6, we obtain two other equations 
for the constants Ci, C2, C3, C4: 

1 + LU 1 — LU 



C16—F 



^,^;i + -;-6)+C26-^F 



-;l-t^;-6 



C3-^Ffi±^,i^;2; ^ 
I — z 



1 — z 



(A36) 



Crd - c.)6^F ( 1 + c.; -6 ] +^2(1+ .)6-'^F f i-^, i-^; 1 - -6 



3(1 -z)2 



C3F(i±^,i^;l ^ 



z - 1 



l-z 



(A37) 



The determinant of the set of homogeneous equations (|A25|) , (|A34|) , ( |A3^ ), and ( |A37| ) gives an equation for 

le can simplify the equations for 
Finally, we obtain the following 



determining the dynamical mass. Since we look for a solution with ^ <C 1, we can simplify the equations for 
Ci by using the corresponding formulas for hypergeometrical functions ^ 
equation for the dynamical mass: 



\/l-LU^A{uj)- 



sinh 


u 




f^i(w)) 


sinh 


UJ 




f 52(^^)) 



3Ci ^ 2 2 ' 



1 fpmexp{h'^) 
log 



(A38) 



where 



A{oj) = 



F{^,^;l+u;; 



_6)F(i-^,i-^;l-^;-6) 



-,1/2 



F{^,'-^-A+u;;-6)F{-^,^-^;l-uj;-6) 



2w ^ F(^^,^-l- 



2 ' 2 
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(A39) 

(A40) 
(A41) 
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and the constant /iq is 



h'^ = 2^(1) - ^ ( ^) ^ f i— 1 ) . (A42) 



In the hmit of large Nf ^ 1 the last equation can be solved explicitly and we find 

™d = ^^exp[-(7V/ + -+7-31og7)], (A43) 

and we used that 1; 2; —6) = log7/6 (note also that i+7 — 31og7w —5.14). It is obvious from comparison 
with Eqs.(|A3|) and ( ]A9[ ) that our solution ( A43| ) satisfies the estimates from below and above, which we obtained 
earlier. 

Notice that up to the preexponential factor the dependence of this solution on Nf coincides with the corre- 



sponding dependence of the (strongly coupling) solution (54) in the case of the reduced QED^+i with a 1-brane. 
The cause of that is the fact that, when Nf ^ 1, in both cases the gauge field propagators are dominated by 
the 1-brane vacuum polarization function, which is the same. All the information about extra dimensions (like 
the number of dimensions, geometry, etc.) is contained in the preexponential factor. 

The reason why the gauge field propagator is dominated by the 1-brane vacuum polarization function in the 



reduced QED2^x rather subtle. An analysis of the gap equation (Al) for the running mass function shows that 
its nontrivial solution is formed by momenta on the interval (m^VG, A): this equation with the low ultraviolet 
cutoff md\/& does not have a nontrivial solution. In the limit Nf ^ 1, the vacuum polarization dominates 
on the interval {md\/6, A), and the equation reduces to a simple Gross-Neveu like equation whose solution is 
TOd exp (-Nf). 

Notice that this is not true for the constant mass approximation, where a nontrivial mass is generated even 
for the low ultraviolet cutoff m^VG- Therefore, the constant mass approximation in this case gives a different 
result for the dynamical mass ( |A9| ) than the correct solution for the running mass function ( A43| ). This is unlike 
the case of the reduced QED3_|_]^ with a 1-brane, where the constant mass approximation is reliable. 



* On leave of absence from Bogolyubov Institute for Theoretical Physics, 03143, Kiev, Ukraine. 
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